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ABSTRACT 

We present a weak-lensing analysis of a region around the galaxy cluster CI 1604+4304 {z = 0.897) on 
the basis of the deep observations with the Hubble Space Telescope {HST) /Wide Field Planetary Camera 
2 (WFPC2). We apply a variant of Schneider's aperture mass technique to the observed WFPC2 field 
and obtain the distribution of weak-lensing signal-to- noise ratio (S/N) within the field. The resulting 
S/N map reveals a clear pronounced peak located about 1'.7 {850h^Q kpc at z = 0.897) southwest of 
the second peak associated with the optical cluster center determined from the dynamical analysis of 
Postman et al. A non-linear finite-field inversion method has been used to reconstruct the projected 
mass distribution from the observed shear field. The reconstructed mass map shows a super-critical 
feature at the location of the S/N peak as well as in the cluster central region. Assuming the redshift 
distribution of field galaxies, we obtain the total mass in the observed field to be 1.0 x lO^^h^f^MQ for 
(z) = 1.0. The estimated mass within a circular aperture of radius 280/i^q kpc centered on the dark 
clump is 2.4 x lO^^h'^^MQ. We have confirmed the existence of the 'dark' mass concentration from 
another deep HST observation with a slightly different (~ 20") pointing. 

Subject headings: cosmology: observations — dark matter — galaxies: clusters: individual (CI 
1604+4304) — gravitational lensing 



1. INTRODUCTION 

Weak shear fields of high-redshift galaxies are promis- 
ing, efficient tools to investigate the mass distribution on 
cluster-supercluster scales (Kaiser & Squires 1993; Lup- 
pino & Kaiser 1997; Kaiser et al. 1998; Bartelmann & 
Schneider 2000) and provide a unique mean to detect dark 
mass concentrations (Schneider 1996; Erben et al. 2000). 
Wide-field weak-lensing surveys of projected mass over- 
densities can probe the statistical clustering properties and 
underlying cosmology (Bahcall et al. 1999). 

Recent observations have revealed the existence of a 
supercluster at a high redshift of z w 0.9 (Lubin et al. 
2000). This supercluster contains two massive galaxy 
clusters, CI 1604+4304 at z = 0.897 and CI 1604+4321 
at z = 0.924. The two clusters are separated by 17' 
on the sky, corresponding to a projected separation of 

9/i^o^ Mpc. C11604+4304 is located at (aj2ooo, '5j2ooo) 
= ( 16h 04m 19.5s, +43° 04' 33'.'9) and one of the optically- 
selected high-redshift cluster candidates studied by Oke, 
Postman, & Lubin (1998). Postman, Lubin, & Oke(1998) 
presented a detailed photometric and spectroscopic survey 
of the cluster and obtained a velocity dispersion of 1226 km 
s~ and a dynamical mass estimate of 6.2 x IO^^/i^q^M©. 

In this Letter, we present a weak lensing analysis of the 
CI 1604+4304 field on the basis of deep images taken with 
the Hubble Space Telescope (HST) /Wide Field Planetary 
Camera 2 (WFPC2). We shall mainly describe the weak- 
lensing signal-to- noise ratio (S/N) analysis of the observed 
field; the details of the mass reconstruction procedure will 
be present elsewhere. Throughout this Letter, we adopt 
f^o = 1 = 0, and Hq = 5O/150 km s"^ Mpc"^; 1' on the 
sky corresponds to 0.52/i^g"^ Mpc at the cluster redshift. 



2. OBSERVATIONS AND DATA REDUCTION 

The cluster CI 1604+4304 field was observed in 1994 and 
1995 using the WFPC2 camera with the F814W filter on 
board the HST (P.I.: J. Westphal; proposal 5234). We re- 
trieved the calibrated data for CI 1604+4304 field from the 
HST archive. Both the 1994 and 1995 observations have 
the same total exposure time of 32 ksec, consisting of 16 
single orbits of 2 ksec. The 1995 pointing covers the central 
region of CI 1604+4304 (Lubin et al. 1998), while the 1994 
pointing is about 20" offset (AY ~ O'.'l, AY ~ 18"2) from 
the 1995 pointing. For this reason, Lubin et al. (1998) dis- 
carded the 1994 data from their analysis. We analyze both 
the 1994 and 1995 data separately. For each observation, 
the data were shifted and combined into the final frame 
to remove cosmic rays using the IRAF/STSDAS task CR- 
REJ. The PC chip was discarded from our analysis because 
of its brighter isophotal limit, so that the final frame for 
each observation consists of three WFC chips. The side 
length of the WFPC2 field is about 2f5 ( l.26h^o Mpc at 
z = 0.897). 

We use the SExtractor package (Bertin & Arnouts 1996) 
for the object detection, photometry, and measurement 
of image shapes. We extract all objects with isopho- 
tal areas larger than 12 pixels (0'.'0996 pixel"^) above 
2(7 pixel" ^ of the local sky level. We calculate the to- 
tal magnitudes /8i4w in STMAG system for each object 
using MAC_BEST parameters in SExtractor. The cor- 
responding detection thresholds for the 1994 and 1995 
data are fisiiw = 26.2 and 26.5 mag arcsec"^, re- 
spectively. We measure the quadruple moments Qij = 

Jd^9d,ejl{9)/ Jd^ei of the surface brightness 1(9) for 
each object, where the center of light is chosen as the 
coordinate origin. From {Qij}, we construct the image 
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ellipticity ^ (Qn - Q22, 2Qi2)/(Tr(Q) + 2(dctQ)i/2) 
(i = 1,2). To construct catalogs of faint galaxies, wc re- 
fer to the following SExtractor output parameters: We 
exclude objects with FLAGS (extraction flags) > 1 from 
our analysis. Wc make star/galaxy separation using the 
CLASS_STAR parameter (stellarity index) and keep the 
objects with CLASS_STAR < 0.2 as faint galaxy candi- 
dates. Objects with FWHMJMAGE (FWHM profile from 
a Gaussian fit to the core) < 4 pixels are excluded since 
image shapes of extremely small objects relative to the 
pixel scale may be affected by the anisotropic PSF. Fi- 
nally, all objects with /8i4w € (25.5, 27.5) are selected as 
background galaxy candidates. This detection and selec- 
tion procedure leads to the final catalogs with total galaxy 
numbers of iVg = 200 and 241, corresponding number den- 
sities of rig = 45.8 and 56.8 arcmin"^, for the 1994 and 
1995 data, respectively. 

3. WEAK LENSING ANALYSIS 
3.1. Strategy 

Our first aim is to obtain the distribution of weak- 
lensing S/N in the data field. Then, the high peaks in the 
S/N maps can be identified as clusters or mass overdensi- 
ties. To do this, we make use of a variant of aperture mass 
statistics. The statistics rely on the fact that the shear 
7 = (71, 72) and the convergence k := S/Scr are related to 
each other through Vk = Dj = where S is the surface 
mass density of the deflector, Scr = {c^ /4:TrG)Ds/ D^D^s 
is the critical surface mass density, and D = {Dij} {i,j = 
1,2) is a differential operator defined by 

fj_[ d/d0, d/d92\ , . 

^-{-0/802 d/d9i) 

(Kaiser 1995). Since u-y is the gradient of k, operating D 
further on yields 

(Ak, 0) = £'^7 = Du-y = {dbfu-y, votu-y). (2) 

In weak lensing limit (k <C 1 and I7I <C 1), the ex- 
pectation value of the image ellipticity, E[e(0)], is the 
shear 7(0). In practice, however, background sources 
have intrinsic ellipticities e^g) , so that weak lensing analy- 
sis involves the smoothing procedure to reduce the noise. 
We denote the smoothed fields by angular brackets ( ): 
e.g., (k) = / (Pe'W{\e -e'\;d)K{e'), where W{e;d) is a 
smooth, continuous window function with a characteristic 
scale of Because of the commutativity between smooth- 
ing and the mass reconstruction (Van Waerbeke 2000), the 
smoothed quantities (/t) and (7) satisfies the same rela- 
tions as those between k and 7: V (k) = D (7) = ■U(-y) ; 
(A(k),0) = £)2^7) = = (div?2(^),rotw(^.)). The 

Laplacian of (k), i.e. divu^^^, is just the convergence 
convolved with the compensated filter function AT4^, and 
hence equivalent to the aperture mass. Defining U(^^^ := 

D (e) , we see that the observable div?!^^^ traces the distri- 
bution of S in the limit of weak lensing, as pointed out by 
Luppino & Kaiser (1997). On the other hand, rotM^^^ mea- 
sures the 'pure' noise. A discretized estimator for Dui-^) is 
given by = (divit^,), rotu^,)) = -n^^^ Y!^=iP{\^ - 



9m\','f))[it{0m]()),^\:{()m]())], whcrc Et and are the tan- 
gential and the radial components of the image elliptic- 
ity e = ei + ie2 defined by 61(0*; 0o) := -3?[e(^)e-2*<^] 
and er(^;^o) := -3[e(^)e^2»0] ^-^^^-^ ^ ^ kxg{e - 0o) and 
p{e;d) = W"{0;i}) - W'{6;^)/e. The noise properties of 
(e) due to the intrinsic source ellipticities are contained 
in the covariant matrix := E[l)2 (e^'^'') -C*^ 

S m=l 

where is the dispersion of the intrinsic source ellipticities 

and we have assumed E[e^^^9jn)e^^H0n)]ij = 5ij5mnO'1/'^- 
The Kronecker delta 5^ in equation (3) ensures that the 
dispersions of div?i^^j and rottT^^) are the same and that the 
two fields are statistically uncorrelated. The local weak- 
lensing S/N v{6) at position 9 is then defined by 

dwU(,){9) ^ Y.raP^\0-^m\;^)et{em;e) 

v[p) := = , = — . 

"^(^^ VS„P^(I^^-C|;^) 

The resulting formula (eq. [4]) is equivalent to the one de- 
rived by Schneider (1996) using aperture mass statistics. 
In the present Letter, we use a Gaussian window function 

of the form WG{9;'d) = exp{—9-^/'d'^)/TT-&'^, in which case 
p{9;^) = 4(6'/i?)2exp(-02/^2)/^^4 ^ pq{9;^); pg{9;^) 

has its maximum at = ?9 and falls off rapidly at 6 > 
These statistics can be applied to the strong-lensing regime 
(e.g., cluster central region), because the contribution of 
image ellipticities to the aperture mass comes mainly from 
galaxies within an annulus at radius ?? and thus we can 
avoid the strong lensing regime. For each independent 
HST observation, we perform a local S/N analysis using 
equation (4). 

To obtain the projected mass distribution, we perform a 
mass reconstruction to the HST/WFC field. Taking into 
account the high redshift (z = 0.897) of CI 1604+4304 
and small field-of-view (2^5 on a side) of the HST /WFC 
field, we adopt a non-linear finite-field inversion method 
developed by Seitz & Schneider (1997), which takes ac- 
count of the source redshift distribution. Since little 
is known about the redshift distribution of field galax- 
ies, we assume a source redshift distribution of the form 
pAz) = /3z2exp[-(z/zo)'^]/r(3//3)zg (Brainerd, Bland- 
ford, & Smaili 1996), in which case the mean redshift (z) 
is given by {z) = ZqT{A/ (3) /T{2>/ (3). In the present Letter, 
we consider only the case ((2;) ,/3) = (1.0, 1.0). 

Once a smoothed ellipticity field (e) [9) is obtained from 
the observed image ellipticities, a convergence map can be 
obtained through the integral equation 

i^oo{9) - «oc = / dPe' H{9, 9') ■ u^^ {9'), (5) 
Ju 

where quantities with 00-subscript represent the values for 
sources at infinite redshift, Koo is the unknown constant 
which represents the average of Kqo within the data field 
hi, H la the kernel which is the gradient of the scalar field 
that satisfies the Neumann boundary problem (see Seitz & 
Schneider 1996), and u^^ = Djoo- In general, the shear 
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is not direct observable, so that the integral equation (5) 
is non-linear and solved iteratively: A mass reconstruc- 
tion scheme applied to the irregular HST/WFG field is 
outlined in Seitz et al. (1996) and Umetsu, Tada, & Futa- 
mase (1999). 

In order to determine the constant Kqc, we employ a 
dynamical mass estimate of CL 1604+4304 obtained by 
Postman et al. (1998). Postman et al. (1998) esti- 
mated the projected mass inside the circular aperture of 
radius i?5oo = SOOft-^Q*^ kpc at the dynamical center to 
be M500 = 0.842 x IO'^^H'^^Mq, corresponding to the 
mean convergence within the circular aperture of Koo.soo = 
^^soo/'i'-Rioo^cr oc = 0.811. However, owing to the small 
field-of-view of the HST/WFPC2, the observed data fields 
do not cover the whole aperture area even for the 1995 
pointing. The effective aperture area covered by the 1995 
data field is about 60% of the circular aperture area. We 
assume that the mean convergence Koo.off inside the effec- 
tive aperture is the same as the mean convergence K00.500 
inside the circular aperture: Kcc,cir ~ Koc,500 = 0.811. We, 
thus, solve the integral equation (5) iteratively varying the 
value of Koo in such a way that the average Koo.eff of the re- 
constructed convergence within the effective aperture area 
equals Ktxi.soo = 0.811. 

3.2. Analysis of CI 1604+4304 

Since the 1995 pointing covers the central cluster region 
where a dynamical analysis was performed by Postman et 
al. (1998), we first analyzed the 1995 data. We calculated 
the fields with d = 0.'4 {150h^^ kpc) on 24 x 24 grid 

points 9ij (except on the PC chip), from which we ob- 
tained the distribution of the weak-lensing S/N within the 
WFC field. Here we have measured the dispersion from 
the rotation field votui^^){9ij) in the following way: We 
calculated the variance a'^^^ of the noise field rot'U(e)(^jj) 
by averaging over the grid points. Equating (7^^^. with the 
average of cr^{9) over the grid points 9ij, we obtained an 
estimate of cTe = 0.21. 

Figure 1 displays the resulting weak-lensing S/N map 
obtained for the 1995 data. We find two S/N peaks above 
a threshold of J^th = 3. The first maximum has a peak 
height of = 4.5 and its location is about I'.l {850h'^Q kpc) 
southwest of the the optical cluster center (marked with 
cross) determined from the dynamical analysis of Postman 
et al. (1998). The peak hight of the second maximum 
is = 3.6 and located in the vicinity of the dynamical 
cluster center. This offset is about O'.lb (75/i^q^ kpc) and 
comparable to the grid spacing of Ofl {^Oh^^l kpc), so that 
the second peak is consistent with the dynamical cluster 
center and associated with the cluster member galaxies 
(Lubin et al. 1998). In this field, there are a total of 12 
cluster members whose redshifts are spectroscopically con- 
firmed by Keck observations (Lubin et al. 1998). However, 
no cluster members are observed at around the first S/N 
peak (see Fig. 4 of Lubin et al. 1998) because objects 
located on the upper edge of the WF4 chip (upper-right 
WF chip) are not contained within the sampled region of 
spectroscopic observations (Postman et al. 1998). 

Although the detection of the mass concentration can- 
not be confirmed spectroscopically from the available data, 
we have checked the stability and reliability of S/N maps. 



We find that the main features, i.e. the locations and 
heights of the peaks, are insensitive to the choice of the 
weight function or filtering scale. It shoud be noted that 
these main features are stable even if objects with large el- 
lipticities (|e| > 0.4) are removed from our galaxy catalog, 
in which case the peak heights of the first and second max- 
ima are 3.7 and 3.0, respectively. There still remains a pos- 
sibility that the detected S/N peaks might be noise peaks 
due to the intrinsic source ellipticities or the anisotropic 
PSF. Van Waerbeke (2000) has shown from the simulated 
data that the probability distribution of the noise peaks 
in mass maps deviates from Gaussian and the probability 
of finding high noise peaks increases owing to the bound- 
ary effects. Since the field-of-view of the HST/WFG is 
small and we adopt a constant filtering scale, the obtained 
S/N maps may be affected by boundary effects. We can 
investigate such possibilities using independent deep ob- 
servation of this field. In Fig. 2, we show the weak-lensing 
S/N map obtained from the 1994 data. Here we have used 
the same filtering scale and grid spacing as for the 1995 
data. Comparing Fig. 2 with Fig. 1, we see that the first 
maximum in the 1994 S/N map corresponds to that in the 
1995 S/N map. The first maximum in the 1994 S/N map 
has a peak height of = 4.2, whose location coincides 
well with that of the first maximum in the 1995 S/N map 
(aj20oo = 16h 04m 15.7s and (5j2ooo = -1-43° 03'47'.'2). The 
identification of the first S /N maximum from two-different 
pointing observations indicates that the S/N maximum 
does not originate from the anisotropic PSF. Moreover, 
since the first maximum in the 1994 S/N map is located 
at the center of the WF4 chip, the S/N measure at the 
first peak cannot be affected by boundary effects. These 
results confirm the detection of 'dark' mass concentration 
associated with the first S/N maximum. 

We have reconstructed the projected mass distribution 
within the 199^ data field. We calculated the smoothed 
ellipticity (e) (9) on 20 x 20 grid points (except on the PC 
chip) using Wg{9;i^) with i) = 0f35 [UGh^^ kpc). In Fig. 
3, we show the reconstructed Koo-map. In mass reconstruc- 
tion, we implicitly assumed that the dark mass concentra- 
tion has the same redshift z = 0.897 as CI 1604+4304. We 
see from Figs. 1 and 3 that the local maxima in the S/N 
maps correspond to the mass peaks and that the weak- 
lensing S/N traces the projected mass, as expected. The 
Koo-map shows a super-critical feature around the dark 
mass concentration (kqo = 1.19) as well as in the clus- 
ter central region. The projected mass within the 1995 
data field is obtained as 1.0 x IO^^H^qAIq. Prom the 
mass map, the enclosed mass in a 280ft-I^Q"'^ kpc radius cen- 
tered on the dark mass concentration is calculated to be 
2.4 X IO^'^/i^q'^Mq: For the mass-to-light ratio, a lower limit 
is obtained as M/L-q « 5OO/150 in solar units from an up- 
perlimit on the dark clump's surface brightness. 

4. DISCUSSION AND CONCLUSIONS 

We have performed a weak-lensing analysis on the deep 
images of a region around CI 1604+4304 taken with the 
HS'T/WFPC2. We detected two significant maxima in the 
resulting weak-lensing S/N map: the second peak associ- 
ated with the dynamical cluster center and the first peak 
located about 1.'7 south of the dynamical center having no 
optical counter parts. The identification of high S/N peak 
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from two-diffcrcnt pointing observations provides a robust 
evidence of dark mass concentration. Assuming that the 
dark mass clump has the same redshift as CI 1604+4304, 
we reconstructed the projected mass distribution in the 
1995 data field, from which we have estimated the mass 
within a circular aperture of radius 280ft.5Q^ kpc centered 
on the dark clump to be 2.4 x lO^^h^^MQ. The obtained 
results indicate that this dark clump could be a compact, 
massive bound system associated with the supercluster at 
z w 0.9; the obtained mass will bo overestimate if the 
redshift of the dark clump is much less than 0.897. 

Our results demonstrate that the weak-lensing S/N 
statistics are powerful and efficient tools even for high- 
redshift cluster surveys (z ~ 1). For a further confirmation 



of the dark mass concentration, deep wide-field observa- 
tions covering the entire snpcrchister field 20') are re- 
quired. Such observations can probe the mass distribution 
in the supercluster field extending over lOh^Q Mpc. 
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Fig. 1. — Contour plot of the weak-lensing signal-to-noise ratio distribution obtained from 241 galaxy images with /8i4w € (25.5, 27.5) taken 
with the HST/WFPC2 in 1995. The solid and dashed lines indicate the positive and negative ones, respectively. The contours are stepped in 
units of 0.5. The filtering scale (0f4) used for the aperture mass measure is marked and the position of the dynamical cluster center is indicated 
(x). The location of the first maximum is marked with -|-. The peak height of the second maximum is 3.6 and its location coincides well with 
the dynamical center. 



6 




X [arcmln] 



Fig. 2.— The same as Fig. 1 but obtained from 200 galaxy images with /sMW € (25.5, 27.5) taken with the HST/WFPC2 in 1994. The 
same coordinate system as in Fig. 1 is used. 
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Fig. 3. — Contour jjlot of the Koo-distributioii reconstructed from the shear field of the 1995 observation. The contours are stepped in 
units of 0.1. The smoothing scale (0'.35) used to calculate the shear field is marked. The shaded area indicates the circular aperture of radius 
500/i^Q kpc centered on the dynamical cluster center (marked with x). The side length is about 2'.5 (1.26/i^q Mpc at z = 0.897). The mass 
distribution is super-critical around the first maximum (marked with +) and the cluster central region. 



